Theory of Neel and Valence-Bond-Solid Phases on the Kagome Lattice of 

Z n-par at acamit e 
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Recently, neutron scattering data on powder samples of Zn-paratacamite, Zn^Cu^^OHQgC^, 
with small Zn concentration has been interpreted as evidence for valence bond solid and Neel 
ordering [l|] . We study the classical and quantum Heisenberg models on the distorted kagome lattice 
appropriate for Zn-paratacamite at low Zn doping. Our theory naturally leads to the emergence 
of the valence bond solid and collinear magnetic order at zero temperature. Implications of our 
results to the existing experiments are discussed. We also suggest future inelastic neutron and 
X-ray scattering experiments that can test our predictions. 

PACS numbers: 



Introduction.- The hallmark of frustrated magnets 
is the existence of macroscopically degenerate classical 
ground states. It is believed that quantum fluctuations 
about such a highly degenerate manifold may lead to 
unexpected quantum ground states. Proposals for emer- 
gent quantum phases include various quantum spin liq- 
uid and valence bond solid (VBS) phases and there has 
been tremendous progress in a theoretical understanding 
of these phases during the last decade Q. On the ex- 
perimental front, ideal materials with spin- 1/2 moments 
(without orbital degeneracy) on frustrated lattices have 
just become available, providing great opportunities for 
testing old and new theoretical proposals [H, 0, H HI • 

One of the prime examples is a series of recent ex- 
periments [3, H, H, 0, [lij E3 on Zn-paratacamite, 
Zn :c Cu4_ a ;(OH)6Cl2, where Cu 2+ ions carry spin- 1/2 mo- 
ments. Most of the attention has focused on the x = 1 
limit Id, but the present paper will address the x = 
limit |l[ when no Zn is present. An important advantage 
of the x = limit is the absence of stoichiometric disor- 
der, which is a serious complication in the interpretation 
of experiments in the x — 1 limit. 

At x = 1 (herbersmithite), the idealized structure 
without stoichiometric disorder, has the Cu moments re- 
siding only on the layered kagome lattices. Remarkably, 
in the experiment no magnetic ordering has been found 
down to 50 mK even though the Curie- Weiss tempera- 
ture is 0(7 w — — 300 K This has raised the hope 
that the quantum ground state of this system may be a 
quantum spin liquid [12, H, 0|- As mentioned, however, 
the presence of stoichiometric disorder makes the inter- 
pretation of the low temperature data a difficult task. 

The situation is very different near x = 0, where there 
is no intrinsic stoichiometric disorder. The magnetic lat- 
tice of the Cu 2+ spin- 1/2 moments form stacks of alter- 
nating (distorted) kagome and triangular lattices. The 
lattice undergoes a structural change around x ~ 0.33; 
monoclinic (rhombohedral) structure for x < 0.33 (x > 
0.33) [HE]]. ^ s a result , the magnetic lattice for x < 0.33 




FIG. 1: Black (blue) represent the Distorted kagome (dual 
dice) lattice. The numbers, in multiples of 1/8, are the frac- 
tional offsets in the height model on the dual dice lattice which 
characterize the Berry's phase effects. 



can be described as weakly-coupled [l| distorted kagome 
lattices (see Fig. Q]for its structure). 

Our theory is motivated by a recent neutron scatter- 
ing experiment on powder samples of Zn-paratacamite 
for small x [l( . These experiments find two phase transi- 
tions at finite temperature for small x. In the low tem- 
perature phase, the neutron scattering data is consistent 
with collinear (or Neel) magnetic ordering while no mag- 
netic ordering is observed in the intermediate phase. Fur- 
thermore, in both phases a heavy gapped spin-1 mode is 
found as well as evidence for dimerization. Finally, a 
weak distortion in the kagome planes is observed at this 
x which vanishes sharply for x > 0.33. Based on these re- 
sults, the authors of Ref. [lj propose that at small x the 
intermediate phase is a VBS phase which then co-exists 
with magnetic ordering in the low temperature phase. 
It should be noted that identification of the intermediate 
phase as VBS ordered is not consistent with the interpre- 
tation of previous /iSR data@,@| that suggests this phase 
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FIG. 2: The large-iV Sp(iV) mean field phase diagram. In 
the small (large) J'/J regime, the Neel (incommensurate) 
long/short range ordered (LRO/SRO) phases arise. The inset 
represents the Neel ordering pattern. 



is magnetically ordered. As such, further experiments, 
preferably on single crystals, providing more direct mea- 
surements of VBS order in the intermediate phase is nec- 
essary to reconcile these experiments. It is our aim to 
provide useful predictions for such studies. 

In the present work, we study the antiferromagnetic 
Heisenberg model on the distorted kagome lattice shown 
in Fig. [T] We consider two inequivalent exchange in- 
teractions, J > J' , consistent with the distortion seen 
in the experiment, where J corresponds to the shorter 
bond-length (see Fig. [1]). We are mostly interested in the 
zero temperature ground states of this model and their 
connection to Zn-paratacamite at x = and possibly 
for small x. Using the Cu-O-Cu angles identified in the 
experiment, one can utilize the Goodenough-Kanamori 
rule to get J'/J m 0.3 [15j. In the classical Heisenberg 
model, we find that energetics chooses the collinear mag- 
netically ordered state shown in Fig. [5] for J'/J < 0.5 
and there exist highly degenerate classical ground states 
for J'/J > 0.5. The collinear state for J'/J < 0.5 has 
precisely the same magnetic order identified in the low 
temperature phase in Zn-paratacamite for small x 

We then investigate the effect of quantum fluctuations 
and possible quantum paramagnetic phases. We use 
the well-documented method of the Sp(N)-generalized 
Heisenberg model where one can change the magnitude 
of "spin" to control the degree of quantum fluctuations 
TJ, ll6| . The large- N mean field phase diagram of this 
model is presented in Fig. [2] Note that one gets the 
same collinear magnetically ordered state for large "spin" 
for J'/J < 0.5-0.8. Understanding of the nature of the 
quantum paramagnetic phase for small "spin" , however, 
requires careful analysis of the spin Berry's phase and 
quantum fluctuation effect about the saddle point solu- 
tion [lil 17 1 . We show that the resulting quantum param- 
agnetic state is a VBS phase depicted in Fig. [3^l. We call 
this the "pin-wheel" VBS state. It is to be distinguished 
from the "columnar" VBS state (shown in Fig. [3Jd) which 
was suggested as a candidate valence bond order in the 




(a)Pin-wheel state 



(b)Columnar state 
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was suggested as a candidate valence bond order m the 
Ref. 11. It can be shown (l6j that the "pin-wheel" state 
can lower the energy by the resonating moves of dimers 



FIG. 3: Pin-wheel and columnar VBS states and their corre- 
sponding lowest energy triplon excitation spectra. 



around the "pin-wheel" structures. 

In order to provide definite predictions for the valence 
bond solid phase, we compute the triplon dispersions 
(shown in Fig. [3fc-d) for both of the VBS phases and 
suggest that inelastic neutron scattering will be able to 
distinguish these phases via their quite different triplon 
dispersions when a single crystal sample becomes avail- 
able. We also suggest that an X-ray scattering experi- 
ment may clearly distinguish the two VBS ordering pat- 
terns via their different further lattice distortions induced 
by the ordering. The expected X-ray structure factors for 
both phases are shown in Fig. 2] 

Classical Heisenberg Model - Possible magnetic order- 
ing patterns in the classical Heisenberg model can be 
investigated by studying the O(N) model in the large- 
N limit [l8[ where the 3 component spin unit vector is 
replaced by an N-component real-valued vector <fi 7 with 
<f>- <p — N. The collinear magnetic order shown in Fig. [2l 
the same magnetic order observed in the experiment, is 
chosen for J'/ J < 0.5 by this method. When J'/J > 0.5, 
a highly degenerate set of wavevectors have the same low- 
est eigenvalue and energetics alone does not determine 
any particular magnetic order. Thus, if magnetic order- 
ing occurs for J'/ J > 0.5 at finite temperatures, it should 
arise via a thermal order by disorder phenomenon. 

Quantum Sp(N) Model and Mean Field Theory - To 
investigate possible magnetically ordered and quantum 
paramagnetic states in the quantum antiferromagnetic 
Heisenberg model, it is useful to generalize the usual spin- 
SU(2) Heisenberg model to an Sp(N) model [lj|- 
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Let us start with the Schwinger boson representation 
of the spin operator S r — b* ra d a p,b r p, where a, (3 =T,|, 
a are Pauli matrices, b ra are canonical boson operators 
and a sum over repeated a indices is assumed. Note 
that we need to impose the constraint — b\ a b ra — 
2S to satisfy the spin commutation relations, where S 
is the spin quantum number. A generalized model is 
then obtained by considering 2N bosonic fields, b ra with 
a = 1,2, ...,2N and the constraint 7i& = b ra b ra = 2SN. 
The simplest such model with Sp(N) symmetry is 



U = - Jw{J a fib\ a b\, ){J l5 b ri b rrs ), 



(1) 



where J a p is a 2N x 2N totally antisymmetric matrix 
which generalizes the Pauli matrix ier 2 from N = 1. 
When k = ni>/N = 2S is fixed while taking the large- 
TV limit, the saddle point solution can be classified us- 
ing both the valence bond order parameter Q rr i — 
{Japb\ a b\., a) and the magnetization induced by a finite 
condensate x ra — (b ra ). The advantage of the large- 
N mean field theory is that we can investigate both the 
large- k semiclassical limit and the small- k extreme quan- 
tum limit on an equal footing (lij. 

In the distorted kagome lattice, we need two inequiv- 
alent valence bond order parameters, Q\ r , and Q^ r ,, as 
depicted in Fig. [T] We also need two Lagrange multi- 
pliers for two inequivalent sites to impose the constraint 
b ra b ra = kN . These parameters need to be determined 
self-consistently in the large-iV mean field theory. The 
large- N Sp(N) mean field phase diagram for the distorted 
kagome lattice is shown in Fig. O For k > k c = 0.26, 
the collinear magnetically ordered state (Fig. [2]) appears 
in the J 1 / J < 0.5-0.8 regime; this is the same magnetic 
order as discovered in the experiment and in the classi- 
cal model. For J'/ J > 0.5-0.8 and at large ft, the ground 
state acquires an incommensurate coplanar order and be- 
comes the \fi x y/3 state at J' = J 2]. The nature of 
the paramagnetic state for small k < n c , however, cannot 
fully be determined within the mean field theory. 

Quantum Fluctuations and Valence Bond Solid in the 
Paramagnetic Phase.- Understanding of the paramag- 
netic phases requires careful analysis of spin Berry's 
phase and quantum fluctuation effects about the mean 
field solution 0, 17 1. It is important to note that 
Q 2 rr , = in the paramagnetic phase for small J'/ J < 0.5- 
0.8. Thus this phase is adiabatically connected to the 
J' = limit, corresponding to the bipartite lattice de- 
picted by the thick lines in Fig. Q] Using Q rr , — Q\e %A ^' , 
one can clearly see that the action in this Q 2 rr , — para- 
magnetic phase is invariant under the U(l) gauge trans- 
formation: b ra e ier b ra (b ra — > e~ l6r b ra ) for r on the 
A-sublattice (B-sublattice) and A rr < — ► A rr < + 9 r — 9 r > ■ 
The effective field theory of such a paramagnetic phase is 
given by the gapped bosonic spinons carrying ±1 gauge 
charges (depending on the sub-lattices) coupled to a U(l) 



gauge field A rr i. Since the spinons are gapped in the 
paramagnetic phase, integrating them out in general pro- 
duces a 2+1 dimensional compact U(l) lattice gauge the- 
ory captured by the simple partition function [2j 

/■ /ll^f exp E v ( cuA p A ) + 1 E ( 2 ) 



where (ij) represent the nearest-neighbor sites of the 
space-time lattice (here we have discretized time) and 
U($) = V(— $) = V($> + 2ir) is an arbitrary periodic po- 
tential. Here p labels the plaquette of the space-time lat- 
tice and curl p A = Yl(ij)<= p s g n P fe)A?' where sgn p (r/) = 
— sgn p = 1 if j comes right after i when one goes 
around a given plaquette p and sgn p (ij) = otherwise. 
Here rjij is an external current determined by spin Berry's 
phase and it is given by rjij = V{rt),{r't') = ±Srr>5t+i,t> 
(for spin-1/2) depending on whether r belongs to the A- 
or B-sublattice. Thus the problem reduces to the com- 
pact U(l) gauge theory with background charges of ±1 at 
the A- and B-sublattice 0]. As well known, this compact 
U(l) gauge theory is confining and the resulting ground 
state would most likely be a VBS. 

In order to find the nature of the VBS state, it is use- 
ful to construct the so-called height model on the dual 
lattice Q, which is equivalent to the compact U(l) gauge 
theory on the direct lattice. The height model can be de- 
rived using the well-documented duality transformation 
and written in terms of the integer-valued height fields 
hj defined on the sites % of the dual space-time lattice 
2\. In our case the dual lattice (in a given time slice) 
is a distorted dice lattice {f} as shown in Fig. Q] (blue 
lattice). Note that the thick blue lines correspond to the 
dual lattice of the J' = limit of the original distorted 
kagome lattice. The height model is found to have action 



(3) 



where (Tj) are the thick bonds of the distorted dual space- 
time lattice and g a non-universal coupling constant. 
Here the constraint hj — hj must be imposed if Tj is 
a thin bond and the offset variables Ci are determined 
by the spin Berry's phase and their time independent 
site-dependent values (i, |, |, 0, |, |, |, |) on the dice 
lattice are shown in Fig. [TJ After solving the simple 
constraint, this height model can be understood using 
standard methods [16|] and the average height fields can 
be determined up to an overall constant [2|. 

The nature of the VBS ground state can be studied 
by using the relation between the height fields and the 
VBS order parameter. It can be shown that the "electric 
field" (in the compact U(l) gauge theory) defined on the 
spatial dual-lattice links is related to the height fields via 
Cfr' = (hf) — (hfi) 0. The VBS order parameter defined 
on the direct-lattice link that intersects the spatial dual- 
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FIG. 4: X-ray structure factor: circles represent Bragg peaks 
of the ideal kagome lattice; triangles arise from the structural 
distortion shown in Fig. [T] These are the only Bragg peaks in 
the columnar state. In the pin-wheel state, additional Bragg 
peaks (hexagons) appear due to further lattice distortion. 

lattice link (ff'} is given by the strength of e f f' 0]- The 
result is a "pin-wheel" pattern shown in Fig. (3^- 

Neutron Scattering: Triplon Dispersion - In the Ref. 
[lj], a different VBS phase (see Fig. [3b) - the "colum- 
nar" phase — was suggested. Here we show that different 
triplon dispersions in the pin-wheel and columnar VBS 
phases can be used to distinguish them if inelastic neu- 
tron scattering experiments are done on single crystals. 

To compute the triplon dispersion, consider letting J 
be the exchange interaction between two spins within the 
same valence bond and A J between two spins on different 
nearby valence bonds. In the decoupled A = limit, the 
triplon dispersion would be completely flat with energy 
J. When A is finite, the triplon band disperses. Here 
we compute this dispersion to first order in A. For this 
purpose, we apply the bond-operator formalism[l!| to 
the valence bonds in the VBS phases where the Hilbert 
space can be represented via singlet and triplet states 
on the bonds of Fig. [3] At first order in A, only the 
processes that preserve triplon number contribute, and 
to this order they become dispersing particles. These 
dispersions for the lowest band in both the pin-wheel and 
columnar states are shown in Fig. [3J:-d. The minima in 
the two cases are clearly located at different positions, a 
feature that can be distinguished experimentally. 

X-ray Scattering- Assuming a lattice contraction 
where valence bonds exist, the pin-wheel state should 
break the lattice translational symmetry of the distorted 
kagome lattice in one of two directions, doubling the unit 
cell. This would lead to new peaks in the X-ray struc- 
ture factor. However, the lattice translational symmetry 
would be intact in the columnar state, leading to no new 
Bragg peaks in the X-ray structure factor in addition to 
those associated with the distorted kagome lattice. The 
X-ray structure factors for the two VBS phases are shown 
in Fig. |U Note that the hexagon symbols represent the 
new Bragg peaks in the pin-wheel state. All other peaks 
also exist in the columnar state. 

Summary and Conclusion - We have provided a theory 
of the zero temperature phases of an antiferromagnetic 
Heisenberg model on a distorted kagome lattice. The re- 
sulting VBS and Neel ordered phases are strikingly sim- 



ilar to those identified in the recent neutron scattering 
experiment on Zn-paratacamite at small doping x [l[ . In 
particular, our theory predicts that "pin-wheel" VBS or- 
dering (see Fig. [3K) can occur as a result of quantum 
disordering of the Neel order. We have suggested future 
neutron and X-ray scattering experiments that can test 
our predictions for this "pin-wheel" VBS ordering. Our 
predictions may also be used for future resolution of the 
disagreement between the interpretations of the neutron 
scattering and /iSR data in the intermediate tempera- 
ture phase. Furthermore, here we focused on zero tem- 
perature ground states of a Heisenberg model so that 
an explanation of the coexistence of Neel and VBS or- 
dering at finite temperature is beyond the scope of this 
work. However, we note that a phase transition between 
the two phases is likely to be first order leaving the pos- 
sibility of a coexist region in the phase diagram. Pos- 
sible relation between the quantum phases on the dis- 
torted kagome lattice described here and the yet-to-be- 
determined quantum ground stated E E Q H ED] 
on the ideal kagome lattice is an important subject of 
future research. 
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SUPPLEMENTARY INFORMATION 

Electromagnetic duality on the distorted kagome lattice. 
The electromagnetic duality of a U(l) lattice gauge the- 
ory on a two dimensional square lattice is well known 
and documented. Here we follow Ref. fsTI ] for the case of 
the duality of a U(l) lattice gauge theory on the strong 
bonds of the distorted kagome lattice of undoped parat- 
acamite. In this 2 + 1 dimensional model, this duality 
mapping connects the U(l) gauge theory with a scalar 
"height" model. An analysis of such a model is then 
straightforward and includes a clear picture for exactly 
how the confinement occurs: the spinons will bind into 
singlets and a valence bond solid (VBS) forms. 
Recall the U(l) lattice gauge theory of Eq. 3: 



TT dA n 



2tt 



exp 



^y(curl p A; • ' X '/<••< 

V (ij) 



(si) 

The first step in transforming to the dual theory is to 
Fourier transform the periodic function exp{V($> p = 
curlp A)} on each plaquette: 

exp V(%) = ex P [V( e p) + ie P® P ] (S2) 



where expV^e p ) is the fourier transform of expV($ p ) 
and e p are integer electric fields living on the plaquettes 
p of the direct lattice. After performing this transforma- 
tion on each plaquette, we integrate out Ay to obtain 

2 = 11 E *(di Vii e- % -W{£ne p )) (S3) 

P {e p £Z} V ' L p J 

where divy is the lattice divergence of e p with plaquettes 
p that share the bond (ij). Thus we have replaced Ay- 
wit h the integer electric fields e p . 

At this stage it is useful to switch to a description on a 
connected dual lattice whose bonds pierce the plaquettes 
of the direct lattice. Such a lattice is actually well known 
in the isotropic limit for the dual of the kagome lattice is 
the dice lattice. Fig. 1 of the main text extends both of 
these lattices to the distorted case. 

The duality mapping is then completed by a suitable 
mapping that relates e p to b-q — — b-p and r/ij to r/ p , where 
1 and j label the sites of the dual lattice of Fig. 1 that 
form the bond (Tj) which pierces the plaquette p of the 
direct lattice and p labels a plaquette of the dual lattice. 
Note: we remove the ambiguity p — ► Tj or p — ► jl by 
associating an outward normal for each plaquette using 
the right hand rule on the A sublattice of the bipartite 



direct lattice. The resulting partition function is 

z = X! IT 5 ( curl p h ~Vp) exp y^Vjhj) 

{{>TjeZ}(T7>6thin P L (TJ)ethick 

(S4) 

Thus the Gauss' law constraint on e p in the direct de- 
scription becomes Ampere's law for bjj in the dual de- 
scription, provided we set bpj = on the thin bonds. 
Aside from this one caveat, this electro-magnetic du ality 
transformation directly follows that on square lattice [Sl|. 

It remains to solve the constraints in Eq. (|S4jl to obtain 
a useful representation of the dual theory. In general, we 
can solve the Ampere's law constraint for bjj by letting 



(S5) 



where h T are integer "height" fields that live on the sites 
of the dual lattice and the integers satisfy 



E(T7> ej5 curl P b ° 



(S6) 



which can take any one of many gauge equivalent con- 
figurations. One possible solution is shown in Fig. ISlh . 
Given such a solution, it is useful to split it up into its 
divergence free and curl free parts: 



(S7) 



where Hjj is the divergence free part and and Hjj are 
fractions whose sum give the integers 6^. For our solution 
Hjj and £i are shown in Fig. ISlb . In conjunction with 
this decomposition it is also useful to specialize to the 
Villain model V(e p ) = — § (e p ) 2 for 

{h T - hj+Q, - Cj+i%) 2 - {h T -hj+b- C 7 ) 2 - e - (Hjjf 

+ terms that vanish in the sum over T ,j (S8) 

so that Hjj only contributes to an overall constant out 
front of the partition function. 

Thus, the solution of the Ampere's law constraint 
leaves us with the height model of Eq. 4 of the main 
text defined on the sites of the dual lattice: 



Z 



IIE 

j h-ez 



n § ( 



(S9) 



with S h = E(rj> e thick H h * ~ h l + C-- Cj) 2 and g a non- 
universal coupling constant. 

Confinement and the pin-wheel VBS phase. The re- 
maining constraints, hj — hj + — if (tj) is a thin 
bond, reduce the number of independent height fields 
and fractional offsets to those living on the sites of a 
square lattice. An understanding of (|S9[) then proceeds 
by a series of mappings onto already known results. 

The dual lattice, shown in Fig. EU has two types of 
sites: black sites with four neighbors which lie on the ver- 
tices of a square lattice and gray sites with two neighbors 
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to* 



(b)docomposition 



FIG. SI: Solution to the Ampere's law constraints on brj on 
the dual lattice (reshaped for convenience), a) one choice of 
the particular solution with fo^ = — = 1 if i and j are 
connected by an arrow from i to j; — otherwise, b) 
The decomposition into divergence and curl free parts: Hjj = 
—H-p = 3/8 if 1 and j are connected by an arrow from 1 to j, 
Hjj = on the dashed line bonds; the fractional offsets are 
the numbers on the sites of the lattice divided by 8. 



by the fractional offsets £„. A simple saddle point anal- 
ysis for this model has already been studied in Ref. 
[S3]. Since ( v takes on four different values depending 
on four different sublattices of the square lattice labeled 
by W, X, Y, Z in Fig. IS21 this saddle point analysis sim- 
ply sets Xv to one of {xw, Xx, Xy , Xz} depending on the 
sublattice and minimizes the energy with respect to these 
four parameters. Following Ref. [S3fl , the result is 



(S13) 
(S14) 
(S15) 
(S16) 
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that are decorations which lie on the bonds of a square 
lattice. The constraint across the thin bonds relates two 
decorations to a vertex. Denoting v as a vertex site and 
d(v) as one of its decorations, we have 



hd(v) — hv, Cd(v) — Cv 



(S10) 



where ( v G {1/8, 3/8, 5/8, 7/8} cyclicly around a plaque- 
tte of the square lattice. Utilizing this mapping, and 
shifting all by 1/8, the height model takes the form 



IlE ex p 



\ - 9 



^ 2 

(vw) 



• y - 

^ 2 

((vw)) 



(h v h w -\- Cw*) 



(Sll) 



with Cv = {0,1/4,2/4,3/4}. This height model differs 
from the simple square lattice case studied by Sachdev 
and Park[Sl| by the addition of off diagonal terms that 
convert it to a model on the anisotropic triangular lattice. 

To gain a simple understanding of the dual height 
model, consider softening the integer height fields, h v — > 
Xv such that Xv is any real number. This necessarily 
introduces periodic potentials, cos27TXu, cos 47rx„, etc. 
The resulting model is then Z x = J T>x e S with 



dr 



E« 2 



(drXv) +2_, (Xv-XwY 
(vw) 



" 2 H (x«-X») 2 + HM 2 cos(2^(x,-C,)) 

((vw)) v 



(S12) 



where for simplicity we have taken the continuum limit 
along the imaginary time axis and we have shifted Xv 



so that the added terms on the diagonal bonds of the 
square lattice do not alter the ground state. 

Given these mean field parameters, we are left with 
the task of determining the physical properties of the 
ground state in terms of quantities defined on the direct 
lattice. To this end, we compute the static magnetic 
fields in the ground state (fey) = Xt — Xj which are the 
gradient of the potential Xt- As stated earlier, each of 
these dual magnetic fields correspond to an electric field 
(hfj) = —(bjj) — ► (e p ) on a plaquette of the direct lattice. 
Spatial plaquettes have a vanishing (e p ), while temporal 
plaquettes have either (e p ) > 0, (e p ) = or (e p ) < 0. A plot 
of (e p ) is shown in Fig. IS2b and has the same symmetry 
as the pin- wheel state shown in Fig. 2 of the main text. 
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(a)4-site 
sublattice 



(b) "Electric" fields 



FIG. S2: (a) The four sublattices W, X, Y, Z of the square 
lattice, (b) the numerator of the square lattice fractional off- 
sets C G {0,1/4,2/4,3/4}. (c) The electric fields (e p ). The 
thickness of the bonds (time-like plaquettes) represent the 
magnitude of (e p ) while the dark bonds have (e p ) > and 
the light bonds have (e p ) < 0. The pattern of dark bonds has 
the same symmetry as the pin-wheel state. 
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